Dynamical chiral symmetry breaking with Minkowski space integral representations 
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The fermion propagator is studied in the whole Minkowski space with the help of the Schwinger- 
Dyson equations. Various integral representations are employed to get solutions for the dynamical 
breaking of chiral symmetry in different regimes of the coupling constant. In particular, in the case 
of massive boson, we extend the singularity structure of the fermion propagator to the two real pole 
Ansatz. 
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In this brief report we investigate a possible scenario 
of dynamical mass generation and estimate the timclikc 
structure of the fermion propagator. This phenomenon, 
dubbed also as dynamical chiral symmetry breaking, re- 
quires intrinsically non-perturbative tools since the parti- 
cle masses can be fully generated via loop contributions. 
In the framework of Schwinger-Dyson equations (SDEs) 
we explore the fermion mass function and the fermion 
propagator in Minkowski space. We develop a novel in- 
tegral technique to solve with reasonable precision a SDE 
which has only been addressed at one-loop order [l[ (for 
Yukawa model). Here, we consider a gauge theory and 
resort to a simple quenched approximation with the mas- 
sive gauge boson transverse mode. The effective coupling 
is then regulated by a Pauli-Villars cutoff A. 

The main result of this paper is to show that for the 
scaling M/A << 1 (walking Technicolor) the analytical 
structure of the exact propagator is given the Lehmann 
representation with one real pole in this propagator. In- 
creasing the ratio M/A, we explicitly show that two pole 
Ansatz plus the corresponding generalized integral rep- 
resentation for the exact propagator is fully adequate for 
the description of dynamical chiral symmetry breaking 
in this phase. The novel integral representation which 
goes beyond the Lehmann representation is introduced 
for this purpose. Within the presented framework we 
achieve larger value of the scaling M/A ~ 0.1. 

In a parity conserving theory the general form of the 
fermion propagator reads 



S(p) 



F( P 2 



p - M(p 2 



(0.1) 



For simplicity we assume F(p 2 ) = 1, which is reasonable 
approximation for gauge theories in the Landau gauge. 
The SDE for the mass function M is modeled in the 
following manner 

M(p 2 ) = S^(p) ~ S-\p) = zg 2 J -0^%)/ 
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where / = p — q, the constant g 2 implicitly absorbs a 
possible group prefactor, and Sq(p) represents the free 



fermion propagator. In such approximation the effective 
coupling does not run logarithmically, but it is constant 
up to a scale A where it rapidly vanishes (i.e. it runs with 
power behaviour). Notice that in what concerns QCD, 
lowering the cutoff to the scale of Aqcd — 250 MeV and 
keeping the coupling large enough such that constituent 
(infrared) quark mass M ~ A can be regarded as an 
approximation of QCD, while when M << A, the limit 
of walking Technicolors is modeled @, lH, for a recent 
treatment within the SDEs framework see 

A little is known about the full Minkowski solution 
of SDEs in strong coupling field theories, hence we can 
refer here the paper of Fukuda and Kugo [6| Furthermore, 
the timelike structure of Greens function as it is read 
from the Euclidean counterparts is not reliably known 
Q. The main aim of this report is to present the direct 
solutions in Minkowski space, assuming a spectral and a 
generalized integral representation of the propagator for 
this purpose. 

In order to carefully compare our Minkowski solutions 
with the spacelike part obtained independently in Eu- 
clidean space, we also perform the Wick rotation iqo — ► q± 
and solve SDE in Euclidean space After the angular in- 
tegration the Euclidean SDE reads [1] , 



M(p%) = 
B(q,p,x) 



3 a 
8^ 



dqj 



M(q 2 E )[B(q,p,m B ) - B(q,p,A)} 



P 2 E l<£ 



M 2 (q E )} 



Ph 



-q 2 E ,x 2 ), 



(0.2) 



where q 2 E = —q 2 (and similarly for p), a = g 2 /(Air) and 
the symbol A stands for the triangle Kallen function, 
A(x, y,z) = (x - y - z) 2 - 4yz . 

The solution of the Eq. (|0.2[) is well known: for the 
coupling below certain critical value a c there exists only 
a trivial solution M(p 2 E ) — 0, while for a > a c we get a 
non-trivial mass function. The value of the critical cou- 
pling a c depends on the details of the kernel, especially 
on the finite ratio R = ms/A, noting that for tub << A 
the critical coupling constant a c — 7r/3 coincides with the 
one obtained in the ladder approximation for the electron 
propagator in the strong coupling QED, where the well- 
known exponential Miransky scaling is exhibited: 



M(0) ~ const ■ Ae 



W(a/a c -l)- 1/2 ] 



(0.3) 
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In the first part of our SDE Minkowski study we as- 
sume spectral representation with a single real pole in the 
propagator and derive the Unitary Equations in their full 
non-linearized form. The solutions of Schwinger-Dyson 
equations obtained by the spectral method has been al- 
ready calculated for several models [1, [t| . Stressed that 
in any case, the resulting spacelike parts of Greens func- 
tions under consideration [g, [^] , were in a good agreement 
with the solutions based on the Euclidean formalism. 

Assumed Lehmann spectral representation reads, 



S(p) = I 

JR+ 



dx 



where S has a pole at p 2 



i>o v {x) 4- a s (x) 



(0.4) 



to 2 , i.e., o s (x) — rm5(x — 
m 2 ) + cr^(x), where r represents the residuum. The func- 
tion cr^(x) is a continuous part of the spectral function 
starting to be non-zero from the first branch point. Sub- 
stituting the integral representation (|0.4|) into our gap 
equation written in Minkowski space, 



M(p 2 ) = iZg 2 \ dx 



diq G(p- q y as{x) 



(27T) 



x + it 



(0.5) 



one arrives to the following dispersion relation for the 
mass function M, 



M(p 2 ) 



duj- 



(4tt) 



dxa s (x) [X (u;m 2 B ,x) - X (u>; A 2 , a;)] , 



Xq{uJ ; a, b) = ^6 (u - {V~a - Vbf) . 

The imaginary part of the propagator and the imag- 
inary part of the dynamical mass function are simply 
related, this relation closes the system of the equations 
(10. 6[) employed. In our approximation Z — 1 and it is 
sufficient to consider the S s part of the propagator, 



Ss(p 2 ) 



M(p 2 ) 



P 



M 2 (p 2 ) 

n 2 -Y 2 - T 2 



(p 2 



i Ej- 



E 2 ) 2 - 



4SIE 2 



(p 2 



E 2 ) 2 + 4S 2 J E 2 



(0.7) 



where we use a shorthand notation, E# = 
Re M(p 2 );E/ = Im M{p 2 ) = -np(p 2 ). Compar- 
ing imaginary part of (|0.7| with the imaginary part of 
the propagator lmS s (p 2 ) = — na^p 2 ), we immediately 
get 



a s V) = Ps {p 2 ) 



p 2 + S 2 J + E 2 



(p 2 ~Z 2 R + Z 2 ) 



4E 2 fl E 2 



(0.8) 



which is nonzero for time- like momenta above the thresh- 
old. The derivation of more general "Unitary equations" 




FIG. 1: Absolute value of the mass function at time-like 
momenta with single pole in the propagator. 



which takes into account the wave function renormaliza- 
tion is straightforward (see for instance 0]). 

The dispersive (real) part of the mass functions is given 
by the principal value integral 



(0.6) 



E fl (;p 2 ) = P ■ Jdtu 



p 2 — to 



(0.9) 



The principal value can be avoided by using p as given in 
(|0.6[) . which yields an ordinary regular integral over the 
new kernel, 

E ^(P 2 ) = £ / dxa s (x) [j(p 2 ,x,m 2 B ) - J(p 2 ,x,A 2 )] , 

where J results from the principal value integration of 
the dispersion relation for M , 



J(p 2 ,x,z) = - 
Q(A P )7(A P ) 



7T p — X — Z 



In 



p — x — z - 



(v^+v^) 2 -p 2 



e(A )V^ 



hi 



p 2 — X — Z + •v/Ap 
An 



-x — z + vAo 



1 



ln(16xz) , 



where we have shortly written A p ( ) = A(p 2 (0), x, z). 

The residuum r and the pole mass function m are ob- 
tained evaluating the dispersion relation and its deriva- 
tive. The coupled set of the integral equations above has 
been solved numerically by iterations. 

The Unitary equations provide solutions for time-like 
momenta above the branch point p 2 > (mg + to) 2 , be- 
low which the propagator is real. The results on the 
negative axis of p 2 are easily obtained by a regular inte- 
gration either from (|0.6[) or from the dispersion relation 



3 



l.Oe+OO - 



I.0e-04 - 




FIG. 2: Dynamical mass function for space-like momenta for 
various coupling a — 1.64; 1.7, 2.0. Solid (dashed) lines stand 
for Euclidean (Minkowski) solutions. 



for M (10. 4p . The resulting timelike solution is presented 
in Fig.rj}. The results presented here are calculated with 
A/to b = 30. 

and we use the mass mj as a scale for all dimensionfull 
quantities. The comparison of spectral Minkowski and 
Euclidean solutions is shown in Fig. [5] Thus, solving the 
Unitary Equations and comparing the Minkowski solu- 
tion to the Euclidean one, we find rather nice agreement 
near the critical coupling. However, when the coupling 
becomes larger (say when a exceeds a c , about a 10 %) 
a discrepancy appears, since the employed spectral rep- 
resentation for the propagator, with just one pole, is no 
longer valid. Retrospectively, the previous one loop ana- 
lytical calculations [l| already found evidence for a more 
complex structure in the propagator. Apparently new 
singularities appear in the propagator for a > ax The 
coupling cvt was determined to be o>t = 1-73 ± 0.02 in 
our case. . 

In what follows we continue with our study of the SDE 
in Minkowski space relaxing our assumption on the spec- 
tral representation of the fcrmion propagator. 

The reason for the fail of one pole Ansatz is easy to 
understand: the dynamical mass is an increasing func- 
tion in the regime from p = to the branch point 
p 2 B = (m + ms) 2 ■ When the coupling is large enough, 
then the large enhancement below the branch point main- 
tains and the dynamical mass function necessarily crosses 
the line \fp 2 ' . Thus, for a certain set of parameters the 
propagator develops two real poles which should be taken 
into account in the integral Ansatz. The second real pole 



first appears at the branch point p B = (m + mj) and 
moves down towards the first one as the coupling in- 
creases. To get a better view we draw this scenario in 
in Fig. [3J The lines displayed represent our numerical 
findings. 

The new integral Ansatz, which is consistent with the 
solution of SDE in the regime where the two poles are 



present, reads 
Ss( P 2 ) -- 



p 2 — m 2 , + iO p 2 — ni 2 _ + iO 
1 



,1, > M .„ . (0.11) 



p 2 — m_ + iO J p 2 — u> + iO 



The alternative equivalent representation can be ob- 
tained by replacements m_ — > — m+, o\T — > a~ . 

This Ansatz immediately implies the dispersive rela- 
tion for the dynamical mass function. Taking F = 1 
approximation again and substituting the formula (jO.llD 
into the SDE it leads after the integration over the mo- 
menta to the following result for the function p 8 {p 2 ) 



i- [X (uj; m 2 + ,m 2 B ) - AT (a>; m\, A 2 ) 



X (uj-,m 2 _,m 2 B ) + X (u;: i m 2 _,A 2 )] 
dx a % ^ [AT (w; m 2 _,m 2 B ) — X (ui; A 2 ,m 2 _)] 



m_ — x 
a s 0) 



dx 



[X {LJ;m B ,x)-X (w;A 2 ,x)] (0.12) 



where the first line in (|0. 1 2[) follows from the first two 
terms in (|0.1ip and the second line follows from the third 
term in (|0.11[) . The derivation of (|0.12p is straightforward 
and follows the same lines as in the case of the standard 
Lehmann representation. The integration over the mo- 
mentum is finite and it remains finite even when A is 
sent to infinity, which is a consequence of the momentum 
behavior of our Ansatz (|0. 1 1 [) . 

The Unitary equations are modified since the integral 
representation has changed. To derive them let us com- 
pare the imaginary and the real parts of the propaga- 
tor above the branch point (m# + m + ) 2 (by definition, 
77i_ > m + ). 

On one side we get from the imaginary part of (|0.11[) 



ImS s (p 2 ) 



7R7 S {p 2 ) 



P 



(0.13) 



and the imaginary part of this function, computed with 
the SDE, is still given by Eq. (f0~7|) . This implies that S s 
is real up to the p 2 = (m_ + tub) 2 ■ 

Thus we get for the time-like momenta, such that p 2 > 
(m_ +%) 2 , the following equation, 



(p - 



m 2 _){p 2 



(0.14) 



This means that the momentum space Schwinger-Dyson 
equation turns into two coupled regular equations (]0 . 1 2[) 
and (|0.14p relating the absorptive and dispersive parts of 
the propagator and its inverse. 

The pole masses are necessarily expected below the 
threshold and they are determined by the zeroes of the 
inverse of the propagator i.e., m± = M(toj_) . The real 
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Our Minkowski solution becomes unstable for a > 2 and 
starts to disagree with the spacelike Euclidean results for 
a > 2. To dive more deep into the chiral breaking phase 
and thus achieve enhancement of the infrared mass would 
require a new reanalysis of up to now unknown, possibly 
complex, propagator singularities. 

To conclude, the one-pole case is adequate to solve 
dynamical chiral symmetry breaking near to the critical 
coupling, while the two pole fit is adequate to work with 
moderate couplings, well into the phase transition region 
where the standard spectral representation deviates from 
the correct solution. For the later suitable integral 

representation of the propagator has been proposed. 
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FIG. 3: (Color online) Propagator singularities bellow thresh- 
old shown for two different coupling constant. The dashed 
(solid) line stands for the case when the propagator function 
Ss exhibits two (one) real poles below the branch points. 

part of the mass entering Eq. (|0.14[) is given by the 
principal value integration over p in Eq. (|0.12[) . This 
leads to the following compact regular integral equation 
for £r, 

- ^B.(p 2 )^ =m_ [J(p 2 ,m^,m%)~ J(p 2 ,m 2 ^A 2 ) 
9 

~J{p 2 , m\,m 2 B ) + J(p 2 ,m\, A 2 )] 

+ / dx^^- \J(p 2 ,m 2 _,m 2 B )- J(p 2 ,m 2 _,A 2 ) 
J m_ — x 

- J(p 2 , x, m%) + J(p 2 ,x, A 2 )] . (0.15) 

The second pole appears for couplings stronger than 
ar — 1.73. In the interval of the couplings 1.73 < 
a < 2.0 the two-pole representation of the propagator 
leads to solutions which, for space-like momenta, agree 
rather well with the Euclidean ones. Our solution for 
the coupling a — 2 is added in the Fig. [2] for spacelike 
and it is displayed in the Fig. 3] for timelike momenta. 
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FIG. 4: (Color online) Absolute value of the mass function at 
time-like regime of the momenta in log- log plot. The upper 
two lines represent the real and imaginary parts of the mass 
function M for a = 2, the two lines below coincide with the 
data plotted in Fig. 2. 
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